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ABSTRACT 

Currently there are a variety of scalar field models that attempt to explain the late time 
acceleration of the Universe. This includes the standard canonical and non-canonical 
scalar field models together with recently proposed Galileon scalar field models. One 
can divide all these scalar field models into two broad categories, namely the thawing 
and the tracker classes. In this work we investigate the evidence for these models with 
the presently available observational data using the Bayesian approach. We use the 
Generalized Chaplygin Gas (GCG) parametrization for dark energy equation of state 
(EoS) as it exhibits both the thawer and tracker like behaviours for different values 
of the parameters. Subsequently we use the Snla ( SN) measurements, the recent 
measurements of the Hubble parameter at different redshifts (H(z)), measurements 
of the lookback time at different redshifts (Lookback), measurements of the linear 
growth factor in large scale structure (GR) from various redshift surveys and finally 
the measurement of the anisotropics in the cosmic microwave background radiation by 
WMAP-7 observations. The analysis of data from SN+H(Z)+Lookback does not favour 
either the tracker or thawer classes. Inclusion of data from GR+ WMAP-7 favours the 
thawer class of models if one assumes the dark energy to be a smooth component. But 
if we consider the dark energy perturbation, both tracker and thawer type of models 
are equally favoured by the data. 
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1 INTRODUCTION 

It is now firmly established that the present Universe is go- 
ing through a phase of accelerated expansion. This has been 
confirmed by different observational results including Type- 
la supernovae ( Riess et al. (1998), Perlmutter et al. (1999), 
Tonry et al. (2003), Knop et al. (2003), Riess et al.(2004), 
Amanullah et al. (2010), Suzuki et al. (2011)), cosmic mi- 
crowave background radiation ( Komatsu et al. (2011)) as 
well as the latest surveys of the large scale structue ( Eisen- 
stein et al. (2005)). The challenge is now to explain the possi- 
ble cause for this late time acceleration and there has been a 
wide range of proposals for this. One of the possible sources 
for this accelerated expansion could be an exotic type of 
matter with an equation of state parameter w < —1/3. Such 
sources are collectively called Dark Energy, [see ( Bean et 
al. (2005), Copeland et al. (2006), Li et al. (2011), Padman- 
abhan(2003), Peebles & Ratra (2003), Salmi & Starobin- 
sky (2000)) for review on dark energy models]. 

The simplest example for this dark energy is a cos- 
mological constant (CC) A, but the need for an extreme 
fine tuning and the cosmic coincidence problem make A 



an unattractive choice from a theoretical perspective. The 
scalar field models are a natural alternative. Different scalar 
field models including quintessence ( Ratra & Peebles (1988), 
Caldwell et al. (1998), Liddle & Scherrer (1999), Steinhardt 
et al. (1999)), tachyon ( Abramo & Finelli (2003), Bagla 
et al. (2003), Aguirregabiria k Lazkoz (2004), Copeland 
et al. (2005)), phantom ( Caldwell (2002)) and k-essence 
( Armendariz et al. (2001), Scherrer (2004), Sen (2006), 
Vikman (2005)) as well as Galileons ( Nicolis et al. (2009), 
Ali et al. (2010), Gannouji & Sami (2010), Appleby & Lin- 
der (2011), Linder (2012), De. Felice & Tsujikawa (2010)) 
have been thoroughly investigated in recent time to explain 
the late time acceleration of the Universe. In addition to 
circumvent the problem like cosmic coincidence, the advan- 
tage of using these fields is that they also mimic CC at late 
times. Moreover, as the equation of state (EoS) for these 
scalar fields in general vary with time, it is expected that 
these scalar field models can be distinguished from ACDM 
with observational data from higher redshifts which are ex- 
pected to become available in near future. 

Scalar field models in general, can be divided into two 
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broad classes: the fast roll and slow roll models, also termed 
as tracking and thawing models, respectively ( Caldwell & 
Linder (2005), Linder (2003)). Tracker models with steep po- 
tentials usually mimic the background matter/radiation in 
the early time and remains subdominant. Only at late times, 
the field leaves the tracker regime and starts behaving like a 
component with large negative pressure and thereby behaves 
like dark energy. On the other hand, thawing models are sim- 
ilar to inflaton, where the scalar field is initially frozen due 
to the large Hubble friction in the flat part of the potential 
and hence behaves close to CC. At late times, the scalar field 
slowly thaws from this frozen state and departs from this CC 
behaviour. In both these cases, a certain degree of fine tun- 
ing (whether in the form of the potentials or in the initial 
conditions) is always necessary to get the required late time 
evolution. 

With the high level of precision that exists in current 
observational techniques, the question naturally arises as to 
whether the available observational data indicate any clear 
preference of one of these models over the other. In the 
present exercise, we attempt to answer this question using 
the Bayesian approach. To do so, it is important to look for 
a parametrization for the dark energy behaviour which in- 
terpolates between the tracker and the thawer behaviour. It 
has been however pointed out that general barotropic mod- 
els described by p = p(p) can not fully mimic the scalar 
field models and in particular the tracker type of models ( 
Scherrer (2006), Linder & Scherrer (2009)). However, the 
attempt in this paper is to have a single model that broadly 
mimics the two overall behaviours, namely: 1) the equation 
of state is frozen initially at w = —1 and then at late times 
it slowly thaws out from this frozen state which is similar to 
thawer type of scalar fields and 2) the equation of state is 
that of matter initially and later on starts behaving like cos- 
mological constant which is similar to tracker type of scalar 
fields. Fortunately a widely studied form, known as " Gener- 
alized Chaplygin Gas" (GCG) ( Kamenshchik et.al. (2001), 
Bento et al. (2002), Sen & Scherrer (2005)), can provide such 
a parametrization scheme [ See also ( Aviles & Cervantes- 
Cota (2011))]. This single parametrization for the EoS of 
the dark energy can not only mimic both of these two be- 
haviours, it can also act as dark energy with constant EoS 
(including uj — — 1 which is the same as CC). Using this 
parametrization we show that considering only the homoge- 
nous background cosmology and applying it to analyze ob- 
servational data does not discriminate between thawing and 
tracker type of behaviours. However, if we include the data 
from the growth of matter perturbations as well as those 
from CMBR anisotropics, thawing type models are preferred 
if we do not consider the perturbation in the GCG fluid. But 
including perturbations in the GCG fluid results no prefer- 
ence of thawer over tracker type of behaviour. 



2 GCG 

We start by discussing the equation of state for GCG given 
by, 

-A 

P(gcg) 8 ' \ 1 ' 
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Figure 1. Behaviour of the equation of state for GCG ^^ gcg ^ as 
a function of redshift. The Solid line is for (1 + (3) < whereas 
the dashed line is for (1 + 0) > 0. A s = 0.85 is assumed. 

where A and /3 are two arbitrary constants ( Bento et 
al. (2002)). Assuming Friedmann- Robertson- Walker metric, 
this equation of state implies that the density evolves with 
redshift as 

P (fl . 3 ,=P (9Cfl ,oh + (l-^)(l + -) 3(1+ ' 3) ] T ^- (2) 
where 

A s = A[ P(gc9)0 ]- (1+m , (3) 

with P( gcg)0 being the present value of P (gcg) • The choice of 
the constants, A s and fi uniquely determines the dark energy 
behaviour in the GCG model. 

One can easily show that A s = —w (gcg) (z = 0). As 
our primary goal is to compare the thawing and tracking 
behaviours which primarily arise in the scalar field models, 
we confine ourselves to non-phantom region i.e A s < 1. For 
A s = 1, GCG exactly behaves as CC at all times and in 
this case the parameter f3 is redundant. For A s — 0, GCG 
behaves exactly as a non-relativistic pressureless matter fluid 
at all times. 

For < A s < 1, the parameter f3 determines the 
how the equation of state changes with time. For f3 = — 1, 
the GCG behaves as a fluid with constant equation state 

P (gcg) = ~ A P(gcg)- For (! + P) > °> at earl y times GCG 

behaves as matter while at late times it behaves as a fluid 
with negative equation of state and asymptotically tends to 
a Cosmological Constant like behaviour in the future. This is 
exactly what we need for a tracking behaviour. On the other 
hand, for (1 + f3) < 0, GCG behaves as a cosmological con- 
stant to start with and then as the Universe evolves, it starts 
deviating from this CC behaviour ( Sen & Scherrer (2005)). 
This is similar to the thawing behaviour. Further, in this 
case, GCG always approaches to matter behaviour assymp- 
totically in the future (z = —1). Hence, in this model, ac- 
celerated expansion is a transient feature. All these different 
behaviours are shown in Figure 1. 

To summarise, GCG behaves as a tracking model for 
(1+/3) > whereas it behaves as thawing model for (1+/3) < 
0. For p — —1, GCG behaves as a fluid with a constant 
equation of state 
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Figure 2. la and 2a confidence contours in the f2 m o — A s — space. The contours are plotted with all the observational data mentioned 
in the text assuming the perturbed dark energy. The upper row is tracker type of models while lower one is for thawer type of models. 



It demarcates the thawing from the tracking regime. 

Considering the fact that GCG can effectively repre- 
sent both the thawing as well as the tracking behaviour by 
suitable parametrization, we use this feature to discriminate 
between these two models on the basis of Bayesian analysis 
using the current observational data. 



3 DATA 

We consider the Supernovae Type la observation which is 
a direct probe for the cosmic expansion. These observations 
measure the apparent luminosity of the Supernovae as ob- 
served by an observer on earth. It is related to the luminosity 
distance cLl(z) defined as 



d L (z) 



(1 + z 



r dz ' 

Jo H(z>) 



(5) 



Using this, we evaluate the distance modulus \i (which 
is an observable quantity) as 



li = m- M 



51 ° g ife + 25 > 



where m and M are, respectively, the apparent and abso- 
lute magnitudes of the Supernovae. We consider the latest 
Union2.1 data compilation ( Suzuki et al. (2011)) consisting 
of 580 data points for the observable /x. 

We also use determinations of the cosmic expansion his- 
tory from red-envelope galaxies by Stern et al. ( Stern et 
al. (2010)). They have obtained a high-quality spectra of 



red-envelope galaxies in 24 galaxy clusters in the redshift 
range 0.2 < z < 1.0 using the Keck-LRIS spectrograph. Sub- 
sequently they have complemented these with high-quality, 
archival spectra as obtained by the SPICES and VVDS sur- 
veys. With this, we get 12 measurements of the Hubble pa- 
rameter H(z) at different redshift s. The measurement of the 
Hubble parameter at present (i.e. z = 0) is obtained from 
the HST Key project. 

The observational data of SNe la and Hubble (H(z)), 
mentioned above relies on the measurements of cosmological 
distances at certain redshifts. On the other hand, the look- 
back time to high redshift objects, as proposed by ( Dalai et 
al. (2001)) depends on the measurement of age of the Uni- 
verse at a particular redshift. The test based on the lookback 
time when combined with the distance measurements can be 
very useful to constrain cosmological models. 

Lookback time is defined as the difference between the 
present age of the Universe (to) and its age at redshift z, t(z) 
and is given by ( Cappzziello et al. (2004)) 



t L (z,p)=t o (0)-t(z) 



(6) = 



1 

Ho 



r dz ' 

J (l + z')U(z>, 



0) 



(7) 

(8) 



where represents the set of model parameters, and 
H(z',9) = H(z',0)/H o . 

The observed lookback time of an object at redshift z% 
is defined as 



t L (z h 6) =t° hs -t(zi) -U 



(9) 



where tg bs is the observed age of the Universe. We use the 
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value of 13.6^q'3 billion years for t° Q hs ( MacTavish et al 
(2005)). The incubation time has been denoted by ti nc and it 
is different for different objects. It is taken to be a nuisance 
parameter in our study and we marginalize over this param- 
eter ( Dantes et al. (2007)). The data for lookback time is 
obtained from ( Simon et al. (2005)). 

We further use the measurement of CMB anisotropy by 
WMAP-7 observations ( Komatsu et al. (2011)). Using the 
publicly available code CAMB ( Lewis et al. (2000)), we 
calculate the expected CMB anisotrpy spectrum for GCG 
as a function of its two model parameters A s and f3 and 
the cosmological parameter Q m o- To incorporate the dark 
energy perturbation, we use publicly available PPF Module 
for CAMB ( Fang et al (2008a), Fang et al. (2008b)). In this 
case, GCG fluid is used as a dark energy with equation of 
state given in eqn (4) and the speed of sound c? s — 1 as this 
fluid represents a parametrization of the scalar field models ( 
For detail calculation and relevant equations in PPF module, 
see ( Fang et al. (2008b))). 

The other cosmological parameters are equated to their 
best fit values as derived from the data from WMAP-7. We 
compare our theoretical results with the data and do the 
likelihood analysis using the likelihood code provided by the 
WMAP team ( Komatsu et al. (2011) ). 

The growth of large scale structures is another useful ob- 
servational probe to distinguish between different dark en- 
ergy models. The growth rate of large scale structures is 
calculated from matter density perturbation 5 = 5pm/ pm in 
the linear regime and subhorizon scales >> H 2 ). On 
these scales where one measures the growth, one can safely 
ignore the dark energy perturbations as it only affects the 
fluctuations comparable to Hubble scales. The effect of dark 
energy is only through the background expansion and one 
can use the Newtonian limit for the evolution of the matter 
density fluctuation satisfying the equation: 

5 + 2-8 -47rGpmS = 0. (10) 
a 

The growth factor / is related to the linear density con- 
trast 5 by / = . The compilation of data for growth 
factor measurements from different redshift surveys can be 
found in ( Gupta et al. (2011)) ( See also references therein). 



4 RESULTS 

Using observational data mentioned above we calculate the 
joint likelihood for the GCG model as: 

- 21og£ = -21og(£ sn x Ch x Ci 00 k x C wm a P x C gr ). (11) 

Before calculating the Bayesian Evidence for thawing 
type and tracker type of models, we maximise the likelihood 
function and calculate the allowed region in the parameter 
space of A s , f3 and Q m o for both thawing type and tracking 
type of models. The results are shown in Figure (2) and 
Figure (3). In these figures we show the 3-D contours in the 
A s — f5— Q m o region. For the tracker type of models ( shown in 
Figure (2 and 3)), it is clear that upto 2 — a confidence level, 
the allowed behaviours are extremely close to C.C. This is 
apparent from the bound on A s which varies negligibly from 
A s = 1, A behaviour. The constraint on Q m o parameter is 



also extremely tight, 0.295 < f2 m o < 0.31. Also, in this case, 
all the values of f3 are allowed. But for two specific values of 
ft, the allowed value of A s is maximum. These are ft = — 1 
for which GCG behaves as a DE with constant equation 
of state and for /3 = for which GCG itself behaves as a 
combination of A and matter fluid. For /3 > only A s — 1 
is allowed which is similar to C.C. 

For thawing case ( shown in Figure (2 and 3)), the vari- 
ation of A s parameter is marginally higher than the tracking 
case, but here also the constraint on Q m o is extremely tight. 
Also in this case, all values of f3 are allowed, but as one 
increases /3, the allowed values for A s get decreased. 

To summarise, with all the data, the GCG parametriza- 
tion puts extremely tight constraint on the dark energy be- 
haviour and the cosmological parameter Q m o. 

We now proceed to substantiate this conclusion. In order 
to do this we calculate the Bayesian Evidence for both the 
thawing as well as tracking behaviour. This is defined as ( 
Liddle et al. (2006)) 

E = J C{0)P{6)d6, (12) 

where represents the set of model parameters, C represents 
the Likelihood function and P(6) is the prior probability dis- 
tribution for parameter 0. In our case, we have three param- 
eters A s , f3 and Q m o . As we have no prior knowledge about 
these parameters, we assume uniform prior for these three 
parameters. We get a thawing behaviour for 1 + f3 < and 
tracking for 1 + f3 > 0. For the thawing case we assume that 
f3 is uniformly distributed in the range — 3 < f3 < — 1 while 
for tracking case we assume that f3 is uniformly distributed 
in the range — 1 < f3 < 1. A s fixes the equation of state 
for dark energy at present (A s = —Wde(z = 0)) and we as- 
sume that it is uniformly distributed between 0.8 and 1 for 
both thawing as well as tracking case. We also assume that 
for both the cases, Q m o is uniformly distributed between 0.2 
and 0.4. 

According to Jeffrey's interpretation ( Jeffreys (1998)), 
if Alni^ between 1 and 2.5, it is a significant evidence for a 
model with higher E while the same between 2.5 and 5 is a 
strong to very strong evidence. If AlnE is more than 5, the 
model with higher E is decisively favoured. 

We calculate AlnE assuming thawing (1 + ft < 0) and 
tracking (1 + f3 > 0) as two different models. We get the 
following results: 

In Ethaw — In Etrac = 0.34 with(sn + hub + look) . (13) 

In Ethaw— In Etrac = —0.55 with(sn+hub+look+gr+wmap) 

(14) 

with D.E. perturbation 

In Ethaw — In Etrac = 1-62 with(sn-\- hub + look + gr + wmap) 

(15) 

without D.E. perturbation 

From these results, it is clear that although observational 
data from the smooth background cosmology cannot dis- 
criminate between thawing and tracking behaviour, when 
these are complemented with the information from the inho- 
mogenous universe but assuming uniform dark energy dis- 
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Figure 3. la and 2a confidence contours in the f2 m o — A s — space. The contours are plotted with all the observational data mentioned 
in the text assuming the smooth dark energy. The upper row is tracker type of models while lower one is for thawer type of models. 



tribution there is significant evidence in favour of thawing 
behaviour ( with odds ratio 5:1). However, when we take 
into account the perturbation in the dark energy compo- 
nent while calculating the CMBR spectra, this preference 
over thawing type of models is lost and there is no Bayesian 
evidence to prefere any type of behaviours. 

We should also mention that although we choose a spe- 
cific prior for A s and Q m o as mentioned above, the results 
are only weakly sensitive to these priors. In order to ver- 
ify this we have varied these priors and we found that the 
above results remain almost the same. This makes our result 
sufficiently robust. 



5 DISCUSSION AND CONCLUSIONS 

Several observed features of the Universe depend on the pre- 
cise nature of cosmological expansion. Broadly the evolution 
of Universe with a dark energy can be classified in two cate- 
gories. In one case the DE equation of state is frozen at the 
w = — 1 in the early universe due to large Hubble friction 
but in the late time as the Hubble friction eases out, the 
equation of state slowly thaws out of this frozen state. This 
is so called thawing behaviour and can be modelled by the 
scalar field with slow-roll potentials. In the other case, the 
DE equation equation of state mimics the background mat- 
ter in the early universe and in the late time it starts behav- 
ing like a C.C. This type of behaviour is called tracker type 
and can be mimicked by the scalar fields with steeper poten- 
tials. Here we parametrize these two overall behaviours by 



the GCG type of equation of state and investigate whether 
the observational data prefer any one of them. 

Towards this end we have first used data pertaining to 
the background cosmology. In particular, we have considered 
three types of data: (i) Supernova Type 1A, (ii) Hubble pa- 
rameter evolution with redshift, and (iii) Lookback time of 
Galaxies. For the Supernova Data we have used the Union2.1 
data. The data from Keck-LRIS, SPICE and VVDS has 
been used for the evolution of hubble parameter while for 
the present redshift hubble parameter HST Key project has 
been used. For the look back time calculation, the incuba- 
tion time of the galaxies has been marginalized over as we 
do not have definitive information about it. 

We have then complemented these with the information 
from inhomogeneous universe. For this we have used the re- 
sults of CMB Anisotropy (as given in the WMAP-7 data) 
and the growth of large-scale structures. 

Our results indicate that if we take into account only 
the smooth background cosmology, then the data is insuffi- 
cient to give a clear verdict either way. On the other hand, 
if we complement this with the information from inhomo- 
geneities in the Universe then there seems to be a preference 
for the thawing model if one does not consider the dark en- 
ergy perturbations while calculating the CMBR spectra. In 
terms of Bayesian evidence, this preference for thawing be- 
haviour is significant. However, we find that incorporating 
the dark energy perturbation to calculate the CMBR spec- 
tra, this preferance for thawing type of behaviour is lost. 

We stress that the GCG cannot represent the full scalar 
field dynamics for tracker type of models in particular and 
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this weakens our conclusion. Nevertheless, our analysis is the 
first one to investigate whether observational data prefer a 
thawing or tracking type of behaviour. 
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